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Abstract. Stars on the lower main sequence (F-type through M-type) have sub- 
stantial convective envelopes beneath their stellar photospheres. Convection in these 
regions can couple with rotation to build global-scale structures that may be observable 
by interferometers that can resolve stellar disks. Here I discuss predictions emerging 
from 3D MHD simulations for solar-type stars with the anelastic spherical harmonic 
(ASH) code and how these predictions may be observationally tested. The zonal flow 
of differential rotation is likely the most easily observable signature of dynamics occur- 
ing deep within the stellar interior. Generally, we find that rapidly rotating suns have 
a strong solar-like differential rotation with a prograde equator and retrograde poles 
while slowly spinning suns may have anti-solar rotation profiles with fast poles and 
slow equators. The thermal wind balance accompanying the differential rotation may 
lead to hot and bright poles in the rapid rotators and cooler, darker poles in slow ro- 
tators. The convection and differential rotation build global-scale magnetic structures 
in the bulk of the convection zone, and these wreaths of magnetism may be observable 
near the stellar surfaces. 



1. Introduction 

When stars like our Sun are young, they rotate rapidly and have strong magnetic fields 
at their surfaces. Magnetic activity is a nearly ubiquitous feature of F- to M-type stars 
on the lower main sequence, all of which have convective envelopes just below their 
photospheres. Younger and more rapidly rotating stars are generally m ore active and 



follow the "rotation-activity" relationship (see e.g.. |Pizzolato et ai]|2003h . Stellar mag- 
netic fields are thought to arise from dynamo processes occuring in these stellar convec- 
tion zones, where turbulent plasma motions couple with rotation to build global-scale 
magnetic fields. At present, stellar dynamo theory does not explain the observed cor- 
relation between rotation and magnetic activity. Likewise, our own Sun's eleven year 
magnetic activity and sunspot cycles remains a puzzling mystery: despite intense study, 
solar dynamo models are at present unable to reliably predict even large-scale features 
of the solar cycle. 

The nature of solar and stellar dynamos, and the origin of solar and stellar magnetic 
fields, remains one of the most important unsolved problems in stellar astrophysics. In 
modern times, the tremendous growth of computational resources, coupled with de- 
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tailed pictures of flows and structure within the solar interior from helioseismology, 
has lead to an explosion of dynamo modeling efforts, ranging from sophisticated 2D 
mean-field models to fully 3D simulations that can capture the non-linear dynamics of 
solar convection self-consistently. Indeed it is now possible to model global-scale con 
vection and dynamo action in the Sun with some fidelity (e .g., Brun & Toomrcl2()() 
Brun et all 120041: iBrowning et all 120061: IMiesch et all l200d l2008l : iMiesch & Toomri . 
2009; Ghizar u et al]l201Ct Racine et all 1201 lh and those efforts are being extended to 



other solar-type stars. The progress in solar simulations was helped tremendously by 
detailed observations of the Sun which constrained and challenged the simulations. 
Detailed observations of other stars will likewise be necessary for further progress in 
a general understanding of stellar convection and dynamo action, which will itself aid 
our understanding of the solar interior. 

Here I will discuss recent global-scale 3D magnetohydrodynamic (MH D) simula- 



tions of stellar convecti on with the anelastic spherical harmonic (ASH) code (IClune et al. 



19991 : iBrun et~alll2004h . This code has been used to simulation solar convection and re- 



produces th e observed solar differential rotation profile relatively well (e.g.. IMiesch et al. 



2006, 2008). Building on this success, we have con ducted a series of s imulations for 



solar- type stars rotating more rapidly than the Sun dBrown et all 12008) an d have ex- 
plored the dynamo-gener ated magnetic fields in several of these cases dBrown et al. 



2010LI2011I :Ib rownll201lh . These simulations are beginning to make specific observa- 



tional predictions which optical interferometry may be able to test, and here we will 
begin laying out what those questions are and how they may be answered. 



2. Simulating stellar convection 

One path towards understanding stellar convection is to conduct simulations of the 
plasma motions occuring within the stellar interior. Stellar convection spans a vast 
range of spatial and temporal scales, which lie well beyond the grasp of direct numer- 
ical simulation even on the largest modern supercomputers. Consequently, models of 
stellar convection and dynamo action must make various tradeoffs, either building up 
from the smallest diffusive scales or building down from the global-scales. The later 
global-scale simulations will be our focus here, and it is these simulations that can 
self-consistently capture the interactions between convection and rotation in a strati- 
fied atmosphere to drive global-scale flows of differential rotation. Th ese simulation s 



remain distant from stellar parameters, and a sense of this gap is given in Brown (201 1). 



3. Global-scale signatures of convection 

The patterns of convection arising in simulations of solar-type stars at a variety of ro- 
tation rates are shown in Figure Q] These simulations span from 0.5 to 10 Q. Q , with 
Q. Q the current solar rotation rate. These are labeled GO. 5 through G10 respectively. 
When the rotation rate is slow relative to convective motions (equivalently, the Rossby 
number is large) then rotational constraints are weak and convective patterns are very 
similar in both polar and equatorial regions. This is evident in case G0.5 (Fig. QJ0. As 
the rotation rate increases (Rossby number decreases), significant differences between 
the equator and pole emerge. Convection near the equator aligns with the rotation axis, 
forming "banana-cells" of convection. Near the poles the convection is more isotropic, 
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Figure 1. Patterns of convection in solar-type stars. Shown are radial velocity 
patterns near the stellar surface for stars rotating at (a) 0.5 Q Q , (b) 3 Q Q , (c) 5 Q Q , 
and (d) 10 Q Q . The broad, slow upfiows are shown in light tones while the narrow, 
fast downflows are shown in dark tones. The north pole is visible and the equator is 
denoted by a dashed line. Clear differences are apparent in the polar and equatorial 
regions, and these become more pronounced as the rotation rate increases. At the 
highest rotation rates, convection near the equator can become confined to narrow 
bands in longitude. To emphasize this in case G10 we show the whole sphere from 
(d) in a Mollweide view in (e), with equator at middle and poles at top and bottom. 
These active nests of convection retain their identity for ma ny thousands of d ays and 
propagate at speeds distinct from the stellar rotation rate dBrown et al.ll2008h . 



with narrow downflow lanes surrounding broader upfiows (e.g., case G5 in Fig.Q]:). As 
the rotation rate is increased the horizontal scale of individual convective cells becomes 
smaller. Individual global-scale convective cells will likely be nearly impossible to de- 
tect on main-sequence solar-type stars; indeed their detection has eluded helioseismic 
detection in the solar interior for many years. 

At the highest rotation rates however, surprising patterns of localized convection 
emerge, and these self-organized structures may create strong observational signatures. 
Here flows near the equator may be confined to one or two active ranges of longitude, 
with quiescent streaming flow in between. One such active nest of convection is shown 
in case G10 in Figure [T}/, with Figure QJ showing the entire near-surface layer in a 
global Mollweide view. These active nests of convection are very long lived structures 
that persist for thousands of days (many hundreds of convective turnover times or rota- 
tion periods). They move at their own angular velocity, distinct from either the stellar 
rotation rate or the differential rotation in which they are embedded and at times may 
cover a substantial fracti on of the stellar dis k. These structures have been found in hy- 
drodynamic simulations dBrown et al.ll2008r) and in some situations they survive in the 
presence of magnetism. In these cases they can act to concentrate surface magnetism 
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Figure 2. Global-scale signatures of convection, (a) Angular velocity of differen- 
tial rotation in rapidly-rotating case GIO. The flow is solar-like, with a fast equator 
and slower flows at the poles. The pole to equator contrast is about 10% near the sur- 
face. This zonal flow is substantially in thermal wind balance and the resulting tem- 
perature profile is shown in (b) with cuts at constant radius shown in (c) (solid near 
surface, dashed mid-convection zone, dotted bottom of convection zone), (d) Black- 
body estimate of change in flux between polar and equatorial regions at the stellar 
photosphere, with hot north pole visible. Scale is relative to stellar average flux and 
is clipped. In slowly rotating stars like case GO. 5, the sense of differential rotation 
may be anti solar (<?) with fast poles and slower equator. The resulting thermal wind 
balance is cold at the poles and warm at the equator (/, g) which may result in dark 
poles (h) due to thermal effects rather than surface magnetism. Here the surface flux 
difference is a few percent. 



into localized structures which may share many similarities with active longitudes of 
stellar magnetism. 

In all solar-like stars the coupling of rotation with convection builds global-scale 
flows of differential rotation and meridional circulation. The profile of differential ro- 
tation is shown for case GIO in Figure |2}?. In the Sun and in more rapidly rotating 
solar-type stars the angular velocity of differential rotation is faster near the equator 
and slower near the polar regions. Thus the equator rotates prograde relative to the 
poles, here with a relative contrast of about 10% between the equator and high-latitudes 
(±60°). In contrast, the meridional circulations are quite weak in these stars, decrease 
in amplitu de with faster ro tation, and are generally multi-cellular in both radius and 
latitude dBrown et al.ll2008h . 

The differential rotation is generally in a thermal-wind balance, and this results 
in a latitudinal temperature structure that is hot near the polar regions and cooler at 
mid-latitudes. In the Sun the magnitude of this contrast is p robably only a f ew Kelvin 
and detecting this signature has been very challenging (e.g. jRast et alj |2008). In more 
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rapidly rotating stars the contrast may be much larger: here in case G10 a contrast of a 
few hundred Kelvin appears in the bulk of the convection zone and likely prints through 
to the stellar surface (Fig. [2^, c). A very simplistic blackbody flux estimate indicates 
that this may cause 10% or larger differences in the surface flux between the polar 
regions and the cooler mid-latitudes (Fig.|2]f). More details on therma l- wind ba l ance i n 
simulations of solar-type s t ars can be found inlBrun & Too mre (20 0J), lMieschl (120051) . 
Miesch et a l. (2006, 2008h. lBallot et all (120071) . and lBrown et all (l2008h . 



At slower rotation rates the dynamics may be substantially different than in the 
rapidly rotating suns. In particular, as the rotational constraint weakens, the sense of 
differential rotation can flip and become anti-solar, with fast poles and a slow equator. 
This is shown for case G0.5 in Figured, where the poles rotate nearly three times faster 
than the equator. These slowly spinning suns can remain in approximate thermal wind 
balance as well, though again the sense flips, now with cool poles and warm equators 
(Fig. Elf, s). 

A striking consequence of this is that the poles may be relatively dark regions, 
due to fluid dynamic effects and irrespective of surface magnetic structures there. Here, 
under a simplistic blackbody flux approximation, the pole is a few percent dimmer than 
the warmer equator (Fig IJfr). Polar spots on slowly rotating stars should be carefully 
examined to see whether any are due to non-magnetic effects. We note briefly that the 
designation "slowly spinning" is dependent on spectral type, and brighter, luminous F- 
type stars will be in the "slow" regime (high Ross by number) even if they rotate several 
times faster than our Sun currently does (lAugustson et al.l 1201 lb . Conversely, less- 
luminous K- and M-type stars are likely "rapid rotators" (low Rossby numbers) even 
when spinning more slowly than the Sun. Dynamo solutions exist for the slowly spin- 
ning suns and these retain an anti-solar differential rotation, though again of reduced 
amplitude relative to the hydrodynamic simulations. 



4. Signatures of internal dynamics 

The scaling of differential rotation and thermal-wind for the rapidly rotating simulations 
are shown in Figure [3] Generally, we find that the angular velocity shear between 
equator and high latitudes AQ grows with faster rotation, though not as quickly as the 
rotation rate itself (a power-law of AQ oc Q 03 is overplotted on cases rotating faster 
than 3 Q Q ). Owing to this, the relative angular velocity contrast 

AQ/Q = (Q(equator) - Q(pole))/Q(equator) (1) 



decreases as the rotation rate increases, here scaling as Q" -° 6 dBrown et al.l2008h . When 



magnetic fields are included and dynamo simulations are conducted, the global-scale 
magnetic fields weaken the large scale differential rotation (Fig.|3j asterisks). At present 
we are still sorting out how magnetism and differential rotation couple, and simulations 
right now are suggesting that the answer may chang e in different regimes of parameter 
space (e.g., cases D10 and D10L in Fig.|3]2, and see lBrownll201 ll) . 



The thermal-wind balance leads to progressively larger latitudinal contrasts of 
temperature in more rapidly rotating simulations (Fig. [3j>). In the most rapidly ro- 
tating cases, the hot poles and the cool mid-latitudes can be nearly 300K different in 
temperature. If this latitudinal temperature gradient prints through the vigorous surface 
convection, then simple blackbody flux arguments would suggest that the poles could 
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Figure 3. Signa tures of differential rotation and the thermal wind (based on 
iBrown et al.ll2008l) . (a) Angular velocity shear of differential rotation AD. in lati- 
tude near the stellar surface shown as a function of rotation rate Q. relative to the 
solar rotation rate Q = 2.6/irad/s. Hydrodynamic cases are shown with diamonds 
while dynamos are labeled and shown with asterisks. AO grows with more rapid 
rotation in hydrodynamic cases. Cases labeled a and b sample more turbulent states. 
(b) Thermal-wind signatures. Shown in blue symbols (and right axis) are the largest 
latitudinal temperature contrasts achieved in these simulations. Shown in black sym- 
bols (and left axis) is the maximum of the relative blackbody flux contrast between 
the bright poles and dim mid-latitudes. 



be up to 20% brighter than the cool mid-latitude bands (Fig. O?, black symbols). The 
magnetic fields in the dynamo simulations reduce this signature as they reduce the dif- 
ferential rotation, and in those simulations the signature is about 5-10%. These black- 
body fluxes are far too simple and neglect nearly all aspects of the stellar atmosphere, 
including viewing angle, but provide an order of magnitude initial estimate. 



5. Wreath-building dynamos 

The magnetic fields produced in these dynamo simulations may also have observational 
consequences. Generally, these rapidly rotating suns generate strong magnetic fields in 
the bulk of their convection zones. This is surprising as many solar dynamo theories 
hold that such organized dynamo action can only occur in the tachocline of shear and 
penetration located between the base of the convection zone and the stable radiative 
zone below. These simulations do not include that interface layer yet they still build 
organized fields. 

Magnetic fields from two of the rapidly rotating dynamo cases are shown in Fig- 
ure |4] These magnetic fields are organized in global-scale wreath-like structures, with 
complex topologies. Magnetic cycles are achieved in many of these dynamo simula- 
tions, including both cases D5 and D10L shown here. In some cases (e.g, D5 in Fig. Hi?) 
the wreaths are highly axisymmetric and have opposite polarities in each hemisphere. 
During the magnetic cycle, wreaths of opposite polarity form in each hemisphere and at 
a later time the polarities will have reversed. Roughly 1500 days later the cycle repeats 
and returns the magnetic fields to a state like is shown here. In some cases (e.g, D10L 
in Fig. the wreaths may be much more concentrated in one hemisphere (here the 
southern). This simulation shows cycles of activity and in each successive cycle the 
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Figure 4. Magnetic wreaths in stellar convection zones, (a) Magnetic wreaths in 
cyclic case D5. Two wreaths of opposite polarity (red, positive; blue, negative) form 
above and below the equator. In this simulation the wreaths undergo reversals of 
polarity on roughly a 1500-day timescale. Relic wreaths from the previous cycle are 
visible in the polar caps, (b) Wreaths in more-rapidly rotating case D10L. Here a 
negative polarity axisymmetric wreath dominates the southern hemisphere while the 
northern hemisphere is filled with non-axisymmetric fields. During the next cycle the 
roles reverse and the wreath instead appears in the northern hemisphere, with tangled 
fields in the southern. Rings of relic fields are strongly evident at the southern pole. 
In both images the colortable saturates at +25kG, while the fields may reach peak 
amplitudes of more than +50kG. 



wreaths alternate between the northern and southern hemisphere, though they rarely fill 
both at the same time. Magnetic fields are present in the other hemisphere (here the 
northern) but are less axisymmetric. In both case D5 and D10L there are significant 
rings of opposite polarit y field located at th e polar caps. These are relic fields from the 
preceding activity cycle (|Brown et al.ll201 ll) . 



6. Constraining simulations with interferometric observations 

Modern simulations of convection in solar-type stars are able to self-consistently gen- 
erate global-scale flows of differential rotation and meridional circulation. In the simu- 
lations, the convection and differential rotation drive strong dynamo action in the bulk 
of the stellar convection zone, generating coherent global-scale wreaths of magnetism. 
The zonal flow of differential rotation is likely the most easily observable signature 
of dynamics occuring deep within the stellar interior. Its characterization is thus of 
cruci al importance. Despi t e this, current obs ervations are in significant disagreement 
(e.g.. [Donahu e et al.lll996l : [Barnes et"aT1l2005l) . and there ar e hints of v ery interesting 
relationships between AQ and the X-ray luminosity of stars ( Saarll2009h . 



Simulations predict that the angular velocity contrast in latitude should grow larger 
in solar- type stars as they rotate more rapidly. At the highest rotation rates, the rela- 
tive contrast AQ./Q. ~ 10%, which at ten times the solar rotation rate (case G10, with 
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Table 1. 


A few possible observational targets 










Name 


HD 


Spectral type parallax 


Vmag 


Vrot 


R/Re 




tet Boo 


126660 


F7V 68 


4.1 


32 


1.5 


10 


chi Dra 


170153 


F7V 124 


3.6 


8 


1.1 


4 


86 Her 


161797 


G5IV 119 


3.4 


8 


1.5 


3 


GJ 702 A 


165341A 


KOV 196 


4.2 


16 


0.8 


10 



rotation period of about 3 days) would lead to an equatorial velocity of about 20 km/s 
and a high-latitude (60° latitude) velocity of about 14 km/s (versus 17 km/s for solid 
body rotation). These are clearly challenging observations but may lie within reach 
for current interferometers (e.g., VEGA at CHARA). To inspire such searches, a few 
of the nearby, rapidly-rotating and bright solar-type stars which should be accessible 
to northern-hemisphere interferometers are listed in Table Q] Rotation rates here are 
lower limits, based on Vsini fro m SIMBAD and rough estimates of stellar radii from 



Pasinetti Fracassini et all (I2001I) . 



The differential rotation may be in thermal wind balance and this may lead to lat- 
itudinal gradients of temperature. These may in term lead to gradients of the stellar 
flux equal to a few percent in brightness. Solar-like differential rotation is probably 
accompanied by relatively warm poles and cool mid-latitudes, while anti-solar differ- 
ential rotation would have cool poles and warm equators. These thermal signatures are 
non-magnetic in nature. 

Lastly, the global-scale longitudinal and radial magnetic fields associated with 
these wreaths may appear at or near the stellar surface. Observations with Zeeman- 
Doppler Imaging (ZDI) techniques appear to show ev idence for large- s cale longitudi- 
nal m agnetic fields at or near the stellar surfaces (e.g. JPetit et al.l l2008; Marsde n et all 
1201 ll ). Interferometric spectroscopic observations may be able to spatially resolve these 
structures, providing critical tests for both the simulations and the growing field of ZDI 
observations, which can in principle be applied to more distant objects. Characteriz- 
ing the large-scale poloidal field may also provide indirect indications of large wreaths 
lurking beneath the stellar photospheres. 
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